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The spectrum of two-dimensional (2D) plasma waves in graphene has been recently studied in the 
Dirac fermion model. We take into account the whole dispersion relation for graphene electrons in 
the tight binding approximation and the local field effects in the electrodynamic response. Near the 
wavevectors close to the corners of the hexagon-shaped Brillouin zone we found new low-frequency 
2D plasmon modes with a linear spectrum. These "inter-valley" plasmon modes are related to the 
transitions between the two nearest Dirac cones. 

PACS numbers: 71.10.-w; 71.45.Gm; 73.21.-b; 73.43.Lp 

Graphene, a recently discovered [3, [3| two-dimensional (2D) material consisting of a single layer of carbon atoms, 
has been in the focus of experimental and theoretical research in the past years (see Ref. [3[ and references therein). 
The carbon atoms in graphene form a dense 2D honeycomb lattice, Fig. [T^,, with two atoms per elementary cell. 
The band structure of graphene electrons 0, [H, [|| consists of two bands touching each other at six points k = K», 
i = 1, . . . , 6 at the corners of the Brillouin zone (BZ, Fig. [TJd). In the vicinity of the dispersion surface forms two 
cones with vertexes at K;. In the intrinsic graphene at zero temperature the lower ("hole") band is fully occupied 
while the upper ("electron") band is empty, and the Fermi level goes through the Dirac points. Using the doping or 
applying a gate voltage between the graphene layer and a substrate (in a typical experiment the graphene layer lies 
on a Si/Si02 substrate) one can shift the chemical potential fj, to the electron or to the hole band and vary the density 
of electrons and/or holes. 

Near the Dirac points the graphene quasi-particles have a linear, quasi-relativistic dispersion 



Etf r = lhV\k-Ki\, |k-Ki|o<l. 



(1) 



Here I = +1 and —1 correspond to the electron and hole band respectively, V ~ 10 8 cm/s is the Fermi velocity in 
graphene and a — 2.46A is the lattice constant, Fig.QJi. It is the massless energy dispersion of graphene quasi-particles 
(fTJ) that leads to its amazing physical properties and caused the great interest to this material. 

In this Letter we address the pro blem of plasma oscillations in graphene. The plasma waves in graphene have been 
considered in Refs. @, @, 0, Hoi. Illj. In these publications the spectrum of plasma waves has been calculated in the 
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FIG. 1: (a) The honey-comb lattice of graphene. All points of the first sublattice (black circles) are given by niai + n 2 a2, of 
the second sublattice (open circles) by niai + n 2 a 2 + b. (b) The BZ of graphene. The basis vectors of the reciprocal lattice 
are Gi and G2. Ki, i — 1, ... ,6 are the corners of the BZ (the Dirac points). Dashed lines show the the boundaries of the 
elementary cell both in direct and reciprocal space. In the figure ai = a(l/2, \/3/2), a 2 = a(— 1/2, v3/2), b = a(0, l/y/3), 
Gi = 27ra _1 (l, 1/V3), G 2 = 27ra _1 (l, -1/Vs), Ki = K 4 = 27ra~ 1 (l/3, 1/V3), K 2 = -K 5 = 2t™- 1 (2/3, 0), K 3 = -K 6 = 
27ra- 1 (l/3,-l/V^), where |ai| = |a 2 | = a, |Gi| = |G 2 | = G = 47r/(^3a). 
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FIG. 2: (a) The Dirac cones in the electron BZ and (b) the low-frequency 2D plasmon modes in the plasmon BZ. 



long- wavelength limit qa <C 1 from zeros of the Lindhard dielectric function 12| , 
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which can be obtained within the self-consistent-field approach (13] or, equivalently, in the random phase approxima- 
tion. Here S is the area of the graphene sample, f(E) is the Fermi-Dirac distribution function, — e is the electron 
charge (e > 0), q is the wavevector of an electric field in the 2D plane, q = |q|, g s = g v = 2 are the spin and valley 
degeneracies and k is the dielectric constant of surrounding medium. The wavefunctions |kZ) have been found from 
the Dirac approach, when the system is described (near the Dirac points) by the effective Hamiltonian Hoir — VcaPa, 
where a takes the values x, y, a a are Pauli matrixes, and p a is the momentum operator. The sub-/superscript "Dir" 
in |j2} reminds that the energies and the wavefunctions have been calculated within the Dirac (effective medium) 
approximation. In the limit q<^kp the spectrum of 2D plasmons takes the form 0,0] 
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which coincides with the standard 2D plasmon dispersion aj p (q) = yj2-Kn Q s e 2 q/ (m*K) with the effective mass being 
replaced by m* — | / V^ 2 . Here kp — \fx\/(HV) is the Fermi wavevector and is the equilibrium surface density of 
charge carriers. At a > kF the curve uj(q) enters the region of the inter-band damping and asymptotically tends to 
the line w = V q 0, S, i, M, El 0] ■ 

The results outlined above are based on Eqs. |T]), and are valid in the "long-wavelength" limit, when both the 
plasmon wavevector q and the Fermi wavevector fcp- are small as compared to the reciprocal lattice vector G ~ 1/a. 
Here we study the 2D plasmon spectrum in graphene at the wavevectors q close to the corners of the BZ. The 2D 
plasmons propagate in the same periodic lattice as the 2D electrons and their spectrum w p (q) should be a periodic 
function of q with the same hexagon-shaped BZ. Near the corners of the plasmon BZ q = K 4 one can expect new 
low-frequency plasmon modes. Indeed, each 2D plasmon wavevector q « corresponds to an intervalley transition 
in the electron BZ Kj — > Kj/, for example, q = Ki corresponds to the transition K5 — > Kg, q = K2 - to the 
transition K@ — ► Ki, and so on. At q sa the energy difference E^i — E^i in the denominator in Eq. (|2|) is close to 
zero, which leads to the new, m£ra-band mfer-valley plasmon modes. In this Letter we show that these low-frequency 
plasmon modes have the linear dispersion, 



Wp(q) = F p |q-K i | ! |q-Ki|a<l, 



(4) 



with the group velocity V p > V. Figure^ schematically shows the low-frequency plasmon mode in the BZ: the central 
dark-grey square-root "flower" and the light-grey "flowers" at the corners of the BZ correspond to the conventional 
mira-valley 2D plasmon ([3]) and the inter-valley plasmons (J4j) respectively. For comparison, the Dirac cones at the 
corners of the electronic BZ are shown in Fig. [5^. 

In order to adequately calculate the graphene response at large wave- vectors q ~ 1/a one should go beyond equation 
@ and take into account the local field effects [15|. The electromagnetic response is then described by the matrix 
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dielectric function 
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e GG > (q, w) = 5 G g> j — -prr Xgg> (q, v) , (5) 

XGG ,(q, c) = |X: /^-^^ (k + q, r|e'^ G >|kO(W| e -^ G )-|k + q, I'), (6) 

where G, G' are reciprocal lattice vectors and xgc is the polarizability tensor. The summation over k in (|5j) is 
performed over the whole BZ, |kZ) are the Bloch functions and E^i is the corresponding energy dispersion. The 2D 
plasmon spectrum is determined by zeros of the determinant of eGG'(q,w), 

det||e GG ,(q,w)|| =0. (7) 

From now on we use the tight binding approximation [B[ for the energy and the wavefunctions. Then the energy 
reads E U i = lA\S k \, where A = 2HV/(y/3a) is the full width of one band and 



5 k = 1 + 2e^ fc W 2 C0S (k x a/2). (8) 



The Bloch functions are 



\kl) = -LVe ,ka 



(CMr-a,2)+I#-a-b,2)), (9) 



where N is the number of elementary cells inside the area S, Ck = <Sk/|<Sk|j ip is the normalized atomic wavefunction 
and z is the perpendicular coordinate. Using the wavefunctions (|9]) one can calculate the matrix elements in ([6]), 

(k/| e -^+ G >|k + q, I') = ijlfflq + G|) [CkCk+q + U'e^ +G ^] , (10) 

where we assume that tp depends only on |r| and M(q) = J d 3 r \ip(r, z)| 2 e* qr . The integration in this formula is 
performed over the whole 3D space. 

If q < 1/a, the terms XGG'/I 1 ! + G l ~ ^(Q + G)M(q + G')/|q + G| in © are small for all G and G' except 
G = G' = 0. Then the general 2D plasmon dispersion equation ([7]) is reduced to the one used in @, d, 0, [13 U|j 
det ||e GG '(q, w)|| = £oo(q, ui) = e£H r (q, a->) = 0. If q is close to one of the vectors K,;, nine terms in the matrix 
XGG'/|q + G| give a noticeable contribution to the determinant of the matrix e G G'(q, uj). For example, if q ~ K 5 , 
the corresponding reciprocal lattice vectors are G = 0, Gi and G2; for them |q+G| ~ |Ks| = K. Thus for q 
close Ki the determinant of the infinite matrix e G G' (q, ui) is reduced to the determinant of a 3 x 3 matrix with 
G, G' = {0,G 1; G 2 }. 

In what follows we assume that q = K5 + q, q <C &f <C K and the temperature T = 0. Then one can show that 
the inter-band contribution (/ 7^ I') to the polarizability tensor $Q is negligible and we can use the linear (Dirac) 
approximation for the energy in the vicinity of both cones. Keeping in ^ only the intra-band terms we get the 
following expression for the polarizability tensor: 

XGG'(q^) = ^M 2 (ir)^ p 6k ~ 9k+ l n *7Gfrq>7G'(fc,q), C 11 ) 

AS ^ E k+ - E k+ q, + + anw 

where a = sgn(/x), 0k = 0(|/-t| — E k +), is the Heaviside step function, ?7 G (k,q) = Ck+qCk + e~ 4 ' q+G ) b and 
G, G' = {0, Gi, G2}; the chemical potential [i can be both positive and negative, i.e. our results are valid for both 
electron and hole gases. Then introducing the notations k = Ki + k, the angle Q\ between the vectors q and k and 
the angle 8 2 between the vector q and K 2 = — K 5 , we get 

3.9, „ 2fT ^ f kF l,l f*" M |k + q|-fc 



4ttW J J -(|k + q| - fc) 2 + {ui/V) 2 

where f G {0) = (-e^ 6 + e~ iGb )/V6 for G = {0,G l7 G 2 } and (f\f) = £ G=0jGl)Ga I/g(#)| 2 = 1- 

If \u>\ < Vq, the integrand in (fT2)) has poles on the d\ axis and the functions xgc turn out to be complex. Similar 
to the standard 2D plasmons (JUS 113 E3 this corresponds to the single-particle intra-band absorption. At larger 
frequencies, \u>\ > Vq, the denominator in (TT2)) does not vanish and the functions Xgg' are real. In this region one can 
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therefore expect a weakly damped (at low temperatures) low-frequency plasmon mode. As it will be seen from the 
result below, the frequency of the new plasmon mode is close to Vq. Evaluating the leading term of the asymptotics 
of Xgg' with respect to the small parameter (|oi|/Vg - 1) < 1 we get 



VH/vq-r 



e GG' - Sgg 7 =^==f G (6 2 )f^(6 2 ), (13) 



where 



0=^=-^-^M 2 (K). (14) 

Noticing that the matrix e can be written as e = 1 — const |/)(/|, and using the formula det (1 — const |/)(/|) = 
1 — const (/|/) we finally get the spectrum of the inter- valley plasmon modes in the form (j4|), where 

V P = V{1 + P 2 ). (15) 

As it is seen from (fT4| the factor (3 is small as compared to unity. 

Calculating the eigenvector corresponding to the inter- valley plasmon ([Jl, (TT3"]) we get the following expression for 
the potential of the plasmon mode 

<f>(r,t) = 4 [f (9 2 )e^ r + / Gl (0 2 )e* Kl - r + /g 2 (0 2 )e lK - r ] ; (16) 

the formula for the corresponding density fluctuation has a similar form. One sees that the inter-valley plasmon with 
the wave vector q = K5 + q is described by a plane wave with the (small) wavevector q and the amplitude periodically 
modulated with the wave vectors Ki, K3 and K5. 

In conclusion, we have found the new infra-band infer- valley low-frequency plasmon modes with the linear dispersion 
(|4]) and the group velocity (jT4j) . The appropriate description of these modes requires to take into account the local 
field effects. The predicted modes do not exist in conventional 2D electron systems and are the unique feature of 



graphene. They could be observed using the Electron energy loss spectroscopy |16|. 
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